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Abstract. We prove that perturbing the reaction-diffusion equation 
Ut = Uxx + {'U-+)'' (p > 1), with time-space white noise produces that 
solutions explodes with probability one for every initial datum, opposite 
to the deterministic model where a positive stationary solution exists. 



1. Introduction 

In this paper we study the following parabolic SPDE with additive noise 
(1.1) ut = u^^ + f{u)+cjW{x,t), 

in an interval (0, 1), complemented with homogeneous Dirichlet boundary 
conditions. Here W is a, 2— dimensional Brownian sheet, cr is a positive 
parameter and / is a locally Lipschitz real function. 

We restrict ourselves to one space dimension since for higher dimensions 
the solution to (|1.1|) (if it exists) it is not expected to be a function valued 
process and have to be understood in a distributional sense. But in this case 
there is no natural way to define /(u), see [17J for more on this. 

Semilinear parabolic equations like (jl.lh arises in the phenomenological 
approach to such different phenomena as the diffusion of a fluid in a porous 
medium, transport in a semiconductor, chemical reactions with possibility 
of spatial diffusion, population dynamics, chemotaxis in biological systems, 
etc. In all these cases, due to the phenomenological approximate character 
of the equations, it is of interest to test how the description changes under 
the effect of stochastic perturbation. 

Equation (jl.ip with / globally Lipschitz has been widely studied (see [17^ 
[T9j). in this case global solutions exist with probability one. However, when 
/ is just locally Lipschitz, typically f{s) ~ with p > 1 or /(s) ~ e*, there 
are practically no results on this problem. Using standard approximation 
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arguments one can easily prove the existence of local in time solutions but 
it does not follow from that proof the behavior of the maximal time of 
existence. 

On the other hand, the deterministic case (i.e. cr = 0) is very well under- 
stood. One problem that has drawn the attention to the PDE community 
is the appearance of singularities in finite time, no matter how smooth the 
initial data is. This phenomena is known as blow-up. What happens is that 
solutions go to infinity in finite time, that is, there exists a time T < oo such 
that 



A well known condition on the nonlinear term / that assures this phenomena 
is when / is a nonnegative convex function with 



For a general reference of these facts and much more on blow-up problems, 
see the book [18] and the surveys [HE]. 

For a large class of nonlinearities /, such as the ones mentioned above, 
problem (jl.ip with a = admits a stationary positive solution v and hence, 
since the comparison principle holds for this equation, for every initial datum 
uq < V the solution to (jl.ip is global in time. 

It is well known (see [6l[18]) that the appearance of blow-up persists under 
(small) regular perturbations. On the other hand, regular perturbations of 
(jl.ip with (7 = admit global in time solutions. Summarizing, the existence 
of global in time/blowing up solutions for this problem with o" = is stable 
under small regular perturbations. Hence it is of interest to test how this 
phenomena is affected by stochastic perturbations. 

Surprisingly, the situation changes for cr > 0. We prove that, in this case, 
there is no global in time solution. In fact, for every initial nonnegative 
datum Uq, the solution to (jl.ip blows up with probability one. 

Stochastic partial differential equations with blow-up has been considered 
by C. Mueller in |14j [T5] and C. Mueller and R. Sowers in [16j. In those 
papers, a linear drift with a nonlinear multiplicative noise is considered and 
the explosion is due to this latter term. 

A similar result, but in some sense in the opposite direction, was proved 
by Mao, Marion and Renshaw in p!3|. There, the authors prove for a system 
of ODEs that arise in population dynamics and that have blow-up solutions, 
that perturbing some coefficients of the system with a small Brownian noise, 
global solutions obtained for every initial data. 

In our problem, a common way to interpret the asymptotic behavior of 
u is the following: consider first the deterministic case a = 0. In this case 
there is some kind of competition between the diffusion, which diffuses the 



lim \\u' 
t/T 



oo 



= OO. 
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zero boundary condition to the interior of the domain and the nonhnear 
source f{u) that induces u to grow very fast. 

Again in the deterministic case, it was proved in that for smah initial 
datum uq, u ^ as t ^ +00, while for uq large, there exists a finite time T, 
such that ||u(-,i)||oo /" +00 as t /" T. More precisely, it is proved that for 
every data uq, there exists a critical parameter A* such that if we solve the 
PDE with initial data Xuq, for A < A* the solution converges to uniformly, 
for A > A* the solution blows-up in finite time and for A = A* the solution 
converges uniformly to the unique positive steady state. 

For small noise a <C 1 one could expect a similar behavior. Of course we 
can not expect convergence to the zero solution as t ^ 00 since in this case 
V = is not invariant for (II. ip . but it is reasonable to suspect the existence 
of an invariant measure close to the zero solution of the deterministic PDE 
and convergence to this invariant measure for small initial datum as t ^ 00. 

However, that is not the case. We prove in Section [3] that for every initial 
datum Uq solutions to (jl.ll) blow-up in finite time with probability one. 

Numerical simulations, as well as heuristical arguments, suggest that, 
for small initial data uq, metastability could be taking place in this case. 
Metastability appears here since, while the noise remains relatively small, 
the solution stays in the domain of attraction of the zero solution of the 
deterministic problem. But, as soon as the noise becomes large, the solution 
escapes this domain of attraction and hence the reaction term begins to 
dominate and pushes forward the solution until ultimately explosion cannot 
be prevented by the action of the noise. 

Organization of the paper. The paper is organized as follows. In Sec- 
tion [2] we give the rigorous meaning of (jl.ip and give the references where 
the foundations for the study of this kind of equation were laid. Section [3] 
deals with the proof of the main result of this paper: the explosion of the 
solutions of (jl.ll) . In Section |4] we propose a semidiscrete scheme in order 
to approximate the solutions to (jl.lh . We prove that the numerical approx- 
imations also explode with probability one and that they converge a.s., in 
time intervals where the continuous solution remains bounded. Finally, in 
Section [5] we show some numerical simulations for this equation. 



2. Formulation of the problem 

We begin this section discussing the rigorous meaning of (II. ip . the refer- 
ences for this being [2| \T7 \ I19j. There are two alternatives: the integral 
and the weak formulation as described in [21 [TTl I19j . The last being more 
suitable for our purposes. Both formulations are equivalent as is shown in 
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Let {J^t)t>o,^) be a probability space equipped with a filtration 

{J^t)t>o which is supposed to be right continuous and such that contains 
all the P— null sets of We are given a space-time white noise on M^- x [0, 1] 
defined on (JS)i>o,IP) and uq £ Co([0,l]). 

Assume for a moment that / is globally Lipschitz, multiply (jl.ip by a test 
function f G C^((0, 1)) n Co([0, 1]) and integrate to obtain 



^0 



.1 .1 
u{x,t)ip{x) dx — / U(){x)ip{x) dx = 
10 Jo 

(2.1) / / u{s,x)(pxx{x) dx ds + / / f{u{s,x))ip{x)dxds 
Jo Jo Jo Jo 

+ cr /" /" ip{x)dW{x,s). 
Jo Jo 

Alternatively, the integral formulation of the problem is constructed by 
means of the function G, the fundamental solution of the heat equation for 
the domain (0, 1). 

u{x,t)- Gt{x,y)uo{y)dy = 
Jo 

t f-l r-t f-l 

I Gt-s{x,y)fiu{y,s))dyds + a / Gt-six,y)dW{y,s). 
Jo Jo Jo 

As a solution to (jl.ip we understand an .Ft— adapted process with values 
in Co([0, 1]) that verifies ([M]) for every ip £ C°°((0, 1)) n Co([0, 1]). 

In [21 [19] it is proved that there exists a unique solution to this problem 
and that the integral and weak formulations are equivalent. 

For / locally Lipschitz globally defined solutions do not exist in gen- 
eral. Nevertheless, existence of local in time solutions is proved by stan- 
dard arguments: consider for each n G N the globally Lipschitz function 
fn{x) = /(-n)l(_oo -„] + /(x)l(_„ „) + f{n)lin,+oo) and n*^, the unique 
solution of (jl.ip with / replaced by /„. Let T„ be the first time at which 
||^^"(') lloo reaches the value n. Then (T„)„ is an increasing sequence of stop- 
ping times and we define the maximal existence time of (jl.ip as T := limr„. 
It is easy to see that u^~^^l^t^Xn} — ^"l{t<r,i} a.s. and hence there exist 
the limit u{x,t) = lim'u"(x,t) for t <T which verifies 

"1 H 
u{x,t f\T)tp{x) dx — I uo{x)ip{x) dx = 
'0 Jo 

/•tAT rl rtAT j-l 

(2.2) / / u{s^x)ipxx{x) dx ds + / / f{u{s,x))(p{x)dxds 
Jo Jo Jo Jo 

/■tAT /■! 

+ cr / / (p{x)dW{x,s). 
Jo Jo 
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So we say that u solves (jl.ip up to the explosion time T. We also say 
that u blows up in finite time if P(T < co) > 0. Observe that if T{u)) < oo 
then 

lim ||ii(-,t,u;)||oo = oo. 



3. Explosions 

In this section, we show that equation (II. 1|) blows-up in finite time with 
probability one for every initial datum uq £ Cq{[0, 1]). Hereafter we assume 
that / is a nonnegative convex function, hence locally Lipschitz. Moreover 
we assume that / 1// < oo. 

In order to prove the blow-up of u, we define the function 

^{t) := / (j){x)u{x,t) dx. 
Jo 

Here (^{x) > is the normalized first eigenfunction of the Dirichlet Laplacian 
in (0, 1). That is, </>(x) = ^ sin(7r2;) and hence we can use it as a test function 
in (j2.ip to obtain 



- $(0) = - Ai /" ^{s)ds+ [ I (t){x)f{u{x,s))dxds 
Jo Jo Jo 



a 



ft fi 

(j){x) dW{x, s). 



Jo Jo 

We denote by zq := $(0) = Jq (I){x)uo{x) dx. 
Now, as / is convex, by Jensen's inequality, we get 

{x)f{u{x, s))dx>f(^l (t){x)u{x, s) dx^ = f{^{s)). 







Moreover, since i;^ is a positive function with L^— norm equal to 1, it is easy 
to see that 

B{t) := — [ [ (t>{x)dW{x,s), 



n 



Jo 



is a standard Brownian motion. 

Combining all these facts, we obtain that $ verifies the (one dimensional) 
stochastic diff'erential inequality 

> ( - Xi^{t) + dt + ^adB{t). 

v8 

Define z{t) to be the one-dimensional process that verifies 
dz = {-Xiz + f{z)) dt + adB, 
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with initial condition z(0) = zq. Then, e{t) = <I>(t) — z{t) verifies 



de > { — Aie H e dt. 

\ ^ - z J 

Observe that e verifies a deterministic differential inequality. Hence, as 
e(0) = it is easy to check that e{t) > as long as it is defined. 

Therefore, > z{t) as long as <I> is defined. 

The following lemma proves that z explodes with probability one. 

Lemma 3.1. Let z be the solution of 

(3.1) dz = {-Xiz + f{z))dt + adB, z(0) = 0. 

Then z explodes in finite time with probability one. 

Proof. The proof is just an application of the Feller Test for explosions ([12], 
Chapter 5). Using the same notation as in [12] we obtain the scale function 
for (EH) to be 



p{x) = j exp(^-^ j b{^)d^^ ds 



Here 6(0 = -Ai^ + /(O- 

It is easy to see that, as f°° 1// < oo, 

p{ — 00) = — oo, p( + Oo) < +00, 

and hence the Feller Test implies that, if S is the explosion time of z, we get 



lim z(t) = +00 = 1 
t/s ' ' J 

To prove that ¥{S < +oo) = 1 we have to consider the function 

I"" p{x) - p{y) 



Jo 0^ 



■dy. 



r'^piy) 

The behavior of v at +oo is given by 1// and hence v(+oo) < +oo, which 
implies that 

¥{S < oo) = 1. 

This completes the proof. □ 

These facts all together, imply that there exists a (random) time T = 
T{uj) < oo a.s. such that 

lim ||n(-, t)||oo = oo a.s. 
So we have proved the following Theorem. 
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Theorem 3.2. Let f be a nonnegative, convex function such that 

-<oo. 

Then, for every nonnegative initial datum uq > Q the solution u to (jl.ip 
blows-up in finite (random) time T with 

p""(r < oo) = 1. 



4. Numerical approximations 

In this section we introduce a numerical scheme in order to compute 
solutions to problem (jl.ip . We discretize the space variable with second 
order finite differences in a uniform mesh of size h = 1/n. That is, for 
X := i/n,i = 1,2, ...,n — 1 the process u'^{t,i/n) = Ui{t) is defined as the 
solution of the system of stochastic differential equations 



(4.1) dui = —^{uij^i — 2ui + Ui^ijdt + f{ui) dt H — — dwi, 2 < i < n — 1, 

h y/h 

accompanied with the boundary conditions ui[t) = Un{t) = 0, Ui{0) = 
uo{ih), 1 < i < n. The Brownian motions Wi are obtained by space integra- 
tion of the Brownian sheet in the interval [ih, (i + l)h). 

Equivalently, this can be written as 

dU = {-AU + f{U)) dt + ^ dW, U{0) = C/°. 

Vh 

Where U{t) = {ui{t), . . . ,u„(t)), —A is the discrete laplacian, f{U) in un- 
derstood componentwise (i.e. f{U)i = f{ui)), dW = {dwi, . . . , dwn) and 
iU% = uo{ih). 

With the same techniques of Theorem 13. 2l it can be proved that solutions 
to this system of SDEs explodes in finite time with probability one. 

We extend n"(t,-) to the whole interval [0,1] by linear interpolation in 
the space variable for each t. 

Concerning the explosions of this system of SDEs we have the following 
Theorem 4.1. Let f be a nonnegative, convex function such that 

oo ^ 

-<oo. 

Then, for every nonnegative initial datum > the solution U to (j4.ip 
blows-up infinite (random) timeT^ with 

p^''(r" < oo) = 1. 
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Proof. The proof uses the same technique of that of Theorem 13.21 Since 
^4 is a symmetric positive definite matrix, we have a sequence of positive 
eigenvalues of j4, < A" < • • • < A^. Let cj)^ the eigenvector associated to 
A". It is easy to see that one can tale cp'^ such that (p^ > for every j, and 
we assume that it is normalized such that Yl^=i — 1- Now, consider the 
function 

n 

i=l 

Proceeding as in the proof of Theorem 13.21 we get that verifies 



where i? is a standard Brownian motion and cr^ air/^/S. The rest of the 
proof follows by Lemma 13.11 as in Theorem 13.21 □ 

Now we turn to the problem of convergence of the approximations. In 
[10) convergence of this numerical scheme for globally Lipschitz reactions is 
proved 

Theorem 4.2 (Gyongy, [lOj Theorem 3.1). Assume f is globally Lipschitz 
and uo G C^{[0,1]). Then 

(1) For every p > 1 and for every T > there exists a constant K = 
K{p, T) such that 

K 

sup sup E(|n"(f,x) -u(t,x)|2P) < — . 
o<i<r2:e[o,i] ""^ 

(2) u^{t,x) converges to u(t,x) uniformly in [0,T] x [0,1] almost surely 
as n —> oo. 

Based on this theorem we can prove that even when / is just locally 
Lipschitz, convergence holds but just in (stochastic) time intervals where the 
solution remains bounded. Observe that a better convergence result is not 
expected. Since the explosion times of u and in general are different, then 
||M'^(f, •) — u{t, •)||oo is unbounded in intervals of the form [0,r] with r close 
to the minimum of the explosion times. To state the convergence result we 
define the following stopping times. Let M > and consider Rm '■= inf{i > 
0, \\uit, OIlL-do,!]) > M} and Rlj := m{{t > 0, ||n"(i, ■)\\l^^[o,i]) > M} 

Theorem 4.3. Assume f is a nonnegative convex function with J j < c«. 
Let u be the solution to (jl.ip and u" its numerical approximation given by 
(liTl) . Then 

(1) For every p > 1 and for every T > there exists a constant K = 
K(p, T) such that 

K 

sup sup E(|u"(t,x) - u{t,x)\'^Pls^t<B^^.I^Bl,}) < ^■ 

0<t<TxG[0,l] 
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(2) For every M > 11^" — ti||Loo([o,TAi?A/]x[o,i]) converges to zero almost 
surely as n ^ oo. 

Remark 4.1. Observe that statement (2) does not make assumptions on the 
numerical approximations u". 

Proof. First, we truncate the / to get a globally Lipschitz function, bounded 
and that coincides with the original / for values of s with \s\ < M. i.e. we 
consider 

' f{s) if \s\ < M 
fMis) = { /(M) iis>M 

fi-M) ifs<-M, 



Let w and be the solutions of (jl.ip and (j4.ip with / replaced by /m 
respectively. 

From Theorem 1321 

K 

sup sup E(|'u;"(t,x) -w{t,x)\^P) < — , 
o<t<Txe[o,i] 

From the uniqueness of solutions of (jl.ip and (j4.ip up to the stopping time 
Rm a R^, we have that almos surely, if t < Rm A ii^^ then u{t, x) = w{t, x) 
and vP'it^x) = w"'{t,x), hence 

sup sup E(|u"(i,x) - j;)|^Pl{t<ij^^/\fl" }) = 

0<t<Tx£[0,l] 

sup sup E{\w''{t,x)-w{t,x)\'^Pl{t<RMAR",})< 

0<t<TxG[0,l] 

K 

sup sup E(|t(;"(t,x) < — . 

0<t<Tx6[0,l] 

This proves (1). To prove (2) observe that since — > w almost surely 
and uniformly in [0, T] x [0, 1] we have that for every e > and < t < Rm, 
||^"(*) OIloo < M + e if n is large enough. That means that liminf R^.^ > Rm 
and hence Rm ^ Rm ~^ R-m- That is the reason we can get rid of K^j. So 
we have 

'"^IIl°°{[0,T]x[0,1]) 
w)l{t<RM-iARlj}\\L^{[0,T]x[0,l]) 
^^)l{t<i?M-iAflJjllL°°{[0,T]x[0,l]) 

^)l{t<RM-i}lli-{[o,r]x[o,i]) 

^llL°°([0,TAi?,Af-i]x[0,l])- 

Since M is an arbitrary constant, this proves (2). □ 



= 


lim 

n— >oo 


- 


> 


lim 

n— >oo 






lim 

n— »oo 


IK- 


> 


lim 

n— >oo 






lim 

n— >oo 
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Remark 4.2. In order to compute an approximate solution this discretization 
in not enough, now we need to discretize the time variable but this is much 
simpler since now we are dealing with a SDE instead of a SPDE. The time 
discretization of (|4.ip can be handled as in 

5. Numerical experiments 

In this section we show some numerical simulations of (II. ip . We perform 
all the simulations with the reaction f{u) = (u^)'^, a = 6.36 and initial 
datum no = 0. 




"4=72'5012 "' " ■■ ' ■ " " "t=7£56"68" " " ' " " " 4=72.5076" 



Figure 1. Profiles of a sample solution at different times. 

To perform the simulations we use the numerical scheme introduced in 
Section [U that is we discretize the space variable with second order finite 
differences in a uniform mesh of size h = 0.02 (i.e.: n = 50 nodes). With 
this discretization we obtain a system of SDE that reads 

dui = Toi'^i+i ~ + Ui-i)dt + f{ui) dt H — p dwi, 2 < i < n — 1, 
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bnapsnot 


lime 


u{-,t,uj) 


oo 


i 


i.UUUU 


5.6159 


2 


50.0000 


3.3863 


3 


72.0202 


15.5104 


4 


72.4202 


18.2885 


5 


72.4802 


38.5848 


6 


72.5002 


82.8705 


7 


72.5012 


203.0799 


8 


72.5068 


2.2695 X 10^ 


9 


72.5076 


1.8128 X 10^^ 



Table 1. The maximum of the solution at differen times 



90 p 
80 - 
70 - 
60 - 
50 - 
40 - 
30 - 
20 - 




Figure 2 . The evolution of the maximum of a sample solu- 
tion with initial data uq = 



accompanied with the boundary conditions ui = Un = 0, Mi(0) = uo{ih), 
1 < i < n. The Brownian motions Wi are obtained by space integration of 
the Brownian sheet in the interval [{i — l/2)h, {i + \/2)h). 

To integrate this system we use an adaptive procedure similar to the one 
developed in [5] for the one dimensional case. Here we adapt the time step 
as in that work replacing the value of the solution (which is a real number) 
by the L^— norm of , as is done in [9] for the deterministic case. More 
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Figure 4. The kernel density estimator of the explosion 
time for a = 6.36 and the corresponding box-plot. 



precisely, the totally discrete scheme reads as follows 
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accompanied with the boundary conditions u\ = Un = 0, for every j > 1 
and = uo{ih), 1 < i < n. Here 

r 

to = 0, Tj = . , tj + 1 - tj = Tj, 

and T is the time-discretization parameter. The Brownian motions Wi are 
the ones of the semidiscrete scheme. 

We want to remark that adaptivity in time is essential in this case since 
a fixed time step procedure gives rise to globahy defined approximations. 

Concerning adaptivity in space, it is knwon for the case a = that it is 
not needed to capture the behavior of the maximal existence time. How- 
ever spatial adaptivity is needed to compute accurately the behavior of the 
solution near the forming singularities (see O El El [9] ) . 

In spite that in Theorem l3.2l we prove that solutions to (jl.ip blow up with 
probability one for every a > and every initial data, we want to remark 
that it is not possible to observe that in numerical simulations since for small 
a, the explosion time is exponentially large when the initial datum is small. 

Essentially, in order to blow-up, the solution needs to be greater than the 
positive stationary solution of the deterministic problem (i.e. the solution 
of 

Vxx — ~f{'^)i which is of size 12 when f{v) — [v^)'^) plus the order of 
the noise a. Once the solution is in that range of values, the noise cannot 
prevent the explosion. 

The probability that such an event occurs in a finite fixed time interval 
depends on a and is exponentially small {pa ~ exp (-l/(j2)). Hence, one can 
estimate P{Tcr > e^/^'^ ) ~ exp(exp(— 1/2(t^)). That means that for a small, 
explosions can not be appreciated numerically and hence the importance of 
the theoretical arguments. 

So, to show the explosive behavior we choose to do the simulations with 
a = 6.36 and initial datum uq = 0. We ran the code with a < 5 until time 
t = 1000 and we did not observe explosions but a meta-stable behavior. 

The features of a particular sample path are shown in Figure [H 

Table [1] shows the times at where the solution is drawn and the L°° — 
norm of the solution at that time. 

In Figure [2] we show the evolution of the L°° norm and in Figure [3] is the 
whole picture as a function of x and t of a sample path. 

Finally, Figure [4] shows some statistics: we perform 832 simulations of the 
solution with a = 6.36 to obtain a sample of the explosion time. Actually, 
we stop the simulation when the maximum of the solution reaches the value 
10^^. The kernel density estimator of the data obtained by the simulation 
and the corresponding box-plot are shown. The sample mean is 46.8834 
and the sample standard deviation 43.8857. 
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These statistics suggest that the distribution of the explosion time T„ is 
close to an exponential variable. This is confirmed by the metastable nature 
of the phenomena. The expected behavior of in this case is 

T 

lim 



.0 E(T^) 

where Z is a mean one exponential variable. 
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